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I. INTRODUCTION 



Building a computer that operates coherently at the quantum level may revolutionise the way we carry out com- 
putations. Indeed it is believed that such quantum computers are much more powerful than their classical analogues. 
For instance it seems that factoring can be carried out exponentially faster on a quantum computer than on a classical 
computer. For this reason much work is being devoted to developing physical systems in which computation can be 
carried out at the quantum level. 

A very attractive systems in which to implement quantum information processing are atomic lattices. Indeed the 
method for realising such lattices suggested in Q has been demonstrated in Q, and lattices comprising more than 10 5 
atoms have been realised. A method for carrying out interactions between neutral atoms suggested in Q has been 
demonstrated in Q. This method realises a global two-qubit gate which in a few steps can entangle all the lattice, 
leading to cluster states^. Finally coherent transport of atoms over many lattice spacings has been demonstrated in 
which implies that the global two-qubit gates can be realised between atoms located many lattice spacings away. 

On the other hand atomic lattices are affected with a fundamental difficulty. Namely it is very difficult in these 
systems to address individually each atom in the lattice. Rather one is limited to the global operations mentioned 
above. Thus whereas atomic lattices seem well suited to carry out simulations of translationally invariant physical 
systems 0, it is not as clear how to use them to implement a universal quantum computer. 

Here we address the question of the computational power of atomic lattices. That is, to what extent can a quantum 
computer be efficiently realised using atomic lattices? 

We shall consider a perfect lattice, i.e. a lattice with exactly one atom per site. We shall suppose that the only 
gates which are available are global one-qubit gates and global two-qubit gates. We will suppose that these gates can 
be performed perfectly. We shall also take each atom to have an internal Hilbert space of dimension 2, i.e. a qubit. 
These restrictions strongly limit the operations that can be carried out and the core of our result consists of showing 
how to overcome these constraint. Finally we shall suppose that the initial state of the lattice breaks slightly the 
translational symmetry in a specific way. Namely we shall suppose that all the atoms are initially in the state |0) 
except two specific atoms that are in the state |1). Our main result is to show that in this situation it is possible to 
efficiently simulate a quantum computer. 

We note that experiments so far involving atomic lattices have only used qubits (as we do), but also have only 
implemented a single global two-qubit Hamiltonian which in the notation below is |01)(01|( d ). On the other hand the 
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result we report here requires two different global two-qubit gates and global one-qubit gates. Whether or not the 
Hamiltonian |01)(01|( d ) and arbitrary global one-qubit gates are enough to simulate a quantum computer is an open 
question. The results obtained here may provide an avenue for tackling this problem. We expect they will also be of 
interest in other contexts as they provide a non trivial way of implementing a quantum computer in a system where 
limited sets of gates are realisable. 

The question of the computational power of atomic lattices has recently been studied in a number of works. For 
instance the proposals of jg and [jj are based on the concept of a "marker qubit" which is circulated through the 
lattice. And uses as ingredient imperfections in the lattice. The latter work has been extended in [Tl| to perfect 
lattices and translationaly invariant initial states. The techniques used in these works are very different from the ones 
presented here. Essential differences concern for instance the size of the Hilbert space of each atom, the initial state, 
and the way local gates between logical qubits are implemented using the global gates. 

Finally it may be interesting to note that the present work was motivated by a numerical study of the computational 
power of atomic lattices. In this numerical investigation we allowed all global one-qubit gates and a single global two- 
qubit gate on qubit pairs of distance 1 on a lattice consisting of n qubits on a circle. We considered the action of 
these global gates on the eigenspaces of the cyclic shift operator in the Hilbert space of the states of the n qubits 
corresponding to the eigenvalue 1. (As the shift operator commutes with the global gates, all of its eigenspaces are 
invariant under the action of global gates.) Using the computational algebra system GAP we obtained that for 
n = 3, . . . , 7 qubits, the restriction of the Hamiltonians of global gates to the eigenspace generate the whole unitary 
Lie algebra. That is, at least up to 7 qubits, the above global gates form a universal set of gates on the eigenspace. If 
this holds for every n (as we conjecture), global one-qubit gates and a single global two-qubit gate on a cyclic lattice 
consisting of n qubits can implement arbitrary unitaries on a Hilbert space of size roughly n — (9(log 2 n) qubits. 

A possible reason for the universality we found is that the global two-qubit gate considered has almost as many 
eigenvalues as possible. But this means that in some sense this gate acts "chaotically". Therefore the model is 
probably not very useful in the sense that it does not seem to allow one to define a qubit structure on the eigenspace 
in a natural way. For this reason we turned to the model described above which uses more two-qubit gates, which 
allows a qubit structure to be defined, and which is amenable to analytic treatment. It is this analysis we report here. 



II. GLOBAL TWO-QUBIT GATES 



We begin by giving a precise definition of global two-qubit gates. 

Let D be a subset of an abelian group G where \D\ = n. The Hilbert space of the pure states of the n qubits is 
C 2 . The elements of the standard basis are indexed by the functions D — > {0, 1}. For a function a : D — ► {0, 1} the 
corresponding basis element is denoted by \a). If p € D we also write a p for the value a(p) € {0, 1}. 

For a 2-qubit operation or 2 2 x 2 2 matrix M and a pair of elements p : q S D, M^ p ' q ^ denotes the n-qubit operation 
which acts as M on the pair of qubits at positions p and q: 

m (p,9) = f M (a p ,a,),(6 p ,b,) if a s = b s for every s G D \ {p,q} : 
a - b 1 otherwise, 

or, in the bra-ket notation 

(a\M M \b) = l ( a p a q\ M \ b p b q) if a s = b s for every s <E D \ {p, q}, 
* ' ' ' [0 otherwise, 

We introduce a weight function W : D x D — > K on the pairs of D. This function corresponds to the fact that 
the global qubit gate can act with different strength on different pairs of atoms in the lattice. We could take W to 
be constant, thereby respecting the translation invariance. For a two-qubit matrix M and a vector d € G the global 
operation is the sum of all copies of M acting on pairs of qubits having difference d, weighted by W: 



AfW>= J2 W(p,q)Mb>«\ 
p,qeD 
p-q = d 



A global 2-qubit Hamiltonian is a matrix of the form where H is an Hermitian (i.e., self-adjoint) 2-qubit oper- 
ation and a global 2-qubit gate is an operation of the form exp(— iH^ ) where is a global two-qubit Hamiltonian. 
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III. AN EFFICIENT ENCODING 

The key to our approach is to use a subset P C D of the qubits as logical qubits. All the qubits in P should 
initially be in a known state, for instance all in the state |0). The rest of the qubits will be set initially to the |0) state 
except for two particular qubits r and r' which are set to |1). The qubits in D \ (P U {r, r'}) play a separator role in 
the computation and are always brought back to |0) after each elementary logical gate, whereas the qubits r and r' 
are always brought back to |1) after each logical gate. The qubits r and r' serve as reference points in our method. 
Intuitively, they are used to "locate" logical qubits and help to "extract" local operations at the right places from 
global ones. The subset P and r, r' obey some geometrical constraints which we now describe. 

1. If p G P, s S {r, r'}, q, q' G P U {r, r'} such that q — q' = p — s then q = p and q' = s. That is for every p G P, 
both p — r and p — r' occur exactly once as a difference of a pair of points from P U {r, r'}. 

2. If q, q' 6 P U {r, r'} such that q — q' = r — r' then q = r and q' = r'. That is r — r' occurs exactly once as a 
difference of a pair of points from P U {r, r'}. 

3. For every p £ P there exists no pair q, q' G P U {r, r'} such that p — r + p — r' = q — q' . 
Examples. It is not difficult to find groups G and sets P, {r, r'} that satisfy these constraints. 

• An / dimensional lattice of size m in each direction: G = Z , D = {0, . . . , m — 1}', P = {p = (pi, . . . ,pi) G D \ 
Y?i=iPi = ( mod 6 )} \ r = (1,0, . ..,0), r' = (2,0,.. .,0). Here |P| = \D\/6-l, i.e., roughly every sixth 
element of D belongs to P. 

• A circle of size n = 6k: D = G = Z„, P = {pEG\p = (mod 6)} \ {0}, r = 1, r' = 2. In this example |P| is 
again \D\/6 — 1. 

• An I dimensional lattice of size m = 3j + 1 in each direction. G = l) , D = {0, . . . ,m — 1} , P — {p — 
(pi, ...,pi) eD\2j + l<p 1 < 3j}, r = (0, . . . , 0), r' = (j, 0, . . . , 0). Here |P| = jm 1 - 1 = ^m 1 - 1 « |£>|/3. 

It is not difficult to generalise these examples or combine them in different ways. 

We say that a function a : D — > {0, 1}™ is admissible if a r = a r t = 1 and a p = for every other p G D \ P. Let 
fc = |P|. Then functions {1, ...,&}—> {0, 1} can be identified with the admissible functions in a natural way therefore 
admissible functions can encode k qubits. 

We can now state our main theorem: 

Theorem 1 Assume that for every pair q ^ q' G P U {r, r'}, we have W(q,q') ^ 0. Let w — max{\W(q, q')\ : q ^ 
q' G D} / min{\W(q, q')\ : g^ij'ePU{r,r'}}. Assume further that for every pair p =/= p' G P , the following global 
two-qubit gates can be implemented for any t and S G {0, 1}: 

exp(-tt(|ll)(ll|)^), exp (-*»(|11><11|)(p-^) , (1) 

exp f-t(\lS){0d\ - \0S){1S\)^ P '^ , exp f-ti{\16) {05\ + \05)(15\) {p - p '^ . (2) 

Then on the Hilbert space of the admissible functions \P\-qubit quantum computations can be efficiently simulated 
using global gates of type 0) and 0). Here by efficiency we mean that the complexity of the simulation, measured in 
the number of global two-qubit gates used, is polynomial in w,n and the complexity of the original computation. 

An upper bound on the efficiency, i.e. on the degree of the polynomial, can be obtained from the proofs of theorems 
^ and [BJ This upper bound is probably far from optimal. 

Section is devoted to the proof of this result. However we shall first show that one can achieve the same result 
as stated in Theorem ^ by using global one-qubit gates and fewer global two qubit-gates. 
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IV. USING FEWER GLOBAL TWO-QUBIT GATES 

The gates eq. Q are the global controlled phase gates. It is interesting to note that the global gates eq. (0) can be 
thought of as generating global Controlled-NOT gates. Indeed the two-qubit Hamiltonians appearing in these gates 
are 

a x ® \8)(8\ , a y <g> \5)(6\ . 

which exponentiated for time tt/2 yield 

exp (-^(|1«5)<(W| + \0S)(1S\)) = (-z(|l><0| + |0)(1|) ® \8)(S\) ; 

and similarly for the other gates in eq. J3J, but in different bases. Note however that the interpretation as a C-NOT 
is not valid for the global gate, because the Hamiltonians acting on the different pairs of qubits do not commute. 

Let us now show that if one can realise global one-qubit gates, then the four global two-qubits eq. can all be 
implemented once a single one can be implemented. To see this we will denote a global one-qubit gate as 



global = -Q u (p) 

peD 



where is the unitary transformation that acts as u on the qubit at position p only: 

(a\u^\b) - 

We then have 



(a p \u\bp) if a s = b s for every seB\ {p}, 
otherwise. 



„, global -iH w O ,globalf u giob a i H (d) global^ 

u e u — e 



exp 



W(p,q)u^u^H^u^u^ 

p,q€D 
p-g = d 



(3) 



Using this expression it is easy to see that the four global 2 qubit-gates appearing eq. @ are equivalent if one can 
implement the global one-qubit gates a 9 x lobal and (1/V5 + ia z /V~2) 9lobal . 

Note that the class of global two-qubit operations that can be implemented when a single global two-qubit operation 
and arbitrary global one-qubit operations can be implemented is larger than the class given in eq. © , see |13| ■ 



V. PROOF OF MAIN THEOREM 

Our aim is to show how local gates between two qubits p,p' £ P can be efficiently implemented by sequences 
of global gates and J3J). We will first study how this can be done at the level of Hamiltonians by commuting 
A(P- r ) = (|ll)(ll|)(P- r ), ^(P-''') = (|ll)(ll|)(p- r ') and the Hamiltonians of the global gates in ©. These results on 
commutation of Hamiltonians will then imply the results for the implementation of two qubit gates, ie. for unitary 
operations. 

For a,b : D — > {0, 1} the elementary matrix with zeros at every position except for a, b where the entry is one is 
denoted by E a y. 



E a ,< 



i)(b\. 



For a : D — * {0, 1} and p £ D we denote by del p (a) the function that can be obtained by zeroing the bit of a at 
position p: 

Our first step will be to investigate how the commutations act on the elementary matrix E a . 
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Recall that A is the two-qubit operation |11)(11|. Its matrix is a diagonal matrix with entry one at position 
corresponding to |11) and zero elsewhere: 



A a ,b = (a\A\b) 



1 if a = b = 11, 
otherwise. 



Thus if 7^ d G G then A^ is the diagonal matrix where the element at the position corresponding to a is just the 
sum of the weights of the 1—1 pairs in a having difference d: 

Ai d l = (a\A^\a)= ]T W(p,q)(a\A^\a) = £ W(p,q). 

p,q G D p,q G D 

p — q = d p — q = d 

a p = a q = 1 

As a consequence, if a is a function from D to {0, 1} and q G D with a q = 1, then we have the following formula: 
^ (d) ^a,del 9 (a)] = (a q - d W(q,q-d)+a q+d W( q + d,q))E aM(t{ay (4) 
(Here, in order to simplify notation, we assume that W(l, I') = if I g" D or I' g" D). Indeed, 

\A {d \E aMi{a) ] = A^\a){del q (a)\-\a){del q (a)\A^ 

£ W(p',q')\a)(de\ q (a)\ - £ W(p', q')\a) (de\ q (a)\. 

p',q'£D p',q'eD 
p' — q' = d p' — q' = d 

a p i = ay = 1 del 9 (a) p / = de\ q (a) q , = 1 

All the terms of the second sum appear also in the first one and the possible terms of the first sum missing from the 
second one are W(q,q— d)\a) (de\ q (a)\ (if a q -d — 1) and W(q + d,q)\a)(de\ q (a)\ (if a q +d = 1)- This gives |@J|. 
A consequence of (@J is that if a is admissible, a q = 1, p G P and s G {r, r'} then 

lA (ps) p 1 _ f W (P> S ) E a.de\ q (a) 11 « G {P' S ) and fl P = !. ^ 



a,del 9 (a)J | o otherwise. 

Indeed, g 6 PU {r, r'} by admissibility of a. If a g _( p _ s ) = 1 then q — (p — s) G PU {r, r'} which implies q = p 
and q — (p — s) = s by the first constraint on PU {?', r'}. In addition g> + p — s G" PU {r, r'} and therefore 
a ?+j9-s — and the coefficient given in |0J is W(p,s). Similarly, a q+p - s = 1 is possible if and only if q = s and 
p = q + p — s and in this case, using once more the first constraint on P U {r, r'} the coefficient is again W(p, s) since 
a p-{p-s)W{p,p — (p — s)) = a s W(p, s) — W(p, s). This discussion proves JSJ. 
On the other hand, if a is not admissible but del g (a) is admissible then 

[^ S) ^a,del,(a)] =°> Provided that {q - p + s, q + p - s} H (P U {r, r'}) = 0. (6) 

This follows from and the fact that all the possible positions where del 9 (a) can be 1 fall in the set P U {r, r'} (by 
admissibility of del 9 (a)). 

These results are the basic ingredients for proving: 

Lemma 2 Assume that a is an admissible function with a q = 1 and b = del g (a). Then for every p G P , we have 

\ A iP-r) \ A (p-t') E ii f W (P> r ) W (P> r ') E a,de\ q (a) if q = P and a p = 1, 

L ' 1 ' o,del,(o)JJ I o otherwise. 

If a is not admissible but del g (a) is admissible then for every p G P, 

[A^\[A^'\E aAeUa) }}=0. 

Proof. Assume that a is admissible. Then repeated applications of (JHJ) (first for s — r' and then for s — r) give that 
[A( p ~ r \ [A( p ~ r \E a ^]] can be nonzero only if q G {p, r'} n {p, r} — {p}, that is q = p. On the other hand, again 
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using JSJ twice, it is straightforward to verify that [A^ p r \[A^ p r \ E a , J] = r) r')E a ,y This 

finishes the proof of the first assertion. 

To see the second statement, assume that the commutator is nonzero. Then, by JHJl, {q — p + r', q + p — r'} fl(PU 
{r, r'}) ^ and {q - p + r, q + p - r} n (P U {r, r'}) ^ 0. Assume first that q - p + r', q - p + r G P U {r, r'}. Then, 
using r — r' = (q — p + r) — (q — p + r') and the second constraint on P U {r, r'}, we have r = q — p + r, r' = q — p + r' 
and q = p G P. But then if a is not admissible then del 9 (a) is not admissible either, a contradiction. (The case 
q + p — r', q + p — r G P U {r, r'} can be treated in a similar way). 

Finally, assume that q+p — r' ,q — p + r *E PL) {r, r'} (the remaining case can be treated by a symmetric argument). 
Then (q + p — r') — (q — p + r) = (p — r) + (p — r'), which is impossible by the third property of the configuration 
PU{r,r'}. □ 

We will now use lemmaElto show how commutations of certain global two-qubit operators Bg P ^ with A^ p ~ r ^ and 

A( p ~ r ' yields a local two qubit operator. The operators Bg P 9 ' will be the basic constituents of the Hamiltonians 
of the global gates in ©. We define B s = \1S}{QS\ for S G {0,1}, i.e. B = |10)(00| and B x = |H)(01|. Assume 
that we take an order of the basis where the first 2' p ' basis elements correspond to the admissible functions and 
the rest correspond to the inadmissible functions. The next lemma states that in this order of basis the matrix of 
[A&- r \ [A<3>- T '\Bf~ q) \] for p,q G P is block diagonal where the upper left 2' p l x 2^ block is a scalar multiple of 
the corresponding block of Bg P ' q \ 

Lemma 3 Let p =/= q G P, a, b : D ~ > {0, 1} such that either a or b is admissible. Then 

(a\{A( p - r \[A^ r '\B { r q) ]}\b)=W(p,q)W^^ 
In particular if only one of a and b is admissible then 

(a\[A {p - r \ {A( p - r '\B < f- q) ]}\b) = 0. 

Proof. For every a, b : D — > {0, 1} and for every pair p' ^ q' G D we have (a\Bg P ' q '\b) — 1 if and only if a p > = 1, 

b p i = 0, a q > = b q i = 5, and a s = b s for every s G D\ {p', q'}. Otherwise (a\Bg P ' q = 0. An equivalent formulation 
of this is 

B g ' 9 = ^ E a,del p ,( a y 

a:D^{0,l}, 
o p ' = 1, a q i = 5 

From this equality we infer 

B { r q) = E W{p',q')Bf^ 

p',q' G D, 
p' -q' =p-q 

E E w(p',i')E aAeUa) 

p',q'eD, a:L>^{0,l}, 
p' - q' = P - q a p > = 1, a q > = 5 

E E W (p'>l') E aAeUay 

a:L>^{0,l} p',q'eD 7 

p' -q' =p-q 
a p i = l,a q i = 6 

Using the latter equality, Lemma El and the fact that [A^ p ~ r \ [A^ p ~ r \E a .b\] is always a scalar multiple of E a ,b, we 
obtain 

(a\\A (p - r) \A {p ~ r "> B {p - q) ]]\b) = I W {P> r ) W {P> r ') W (P><l) if b = delp(a), a p = 1 and a q = 6, 
^ ^ ' ^ 1 ' [0 otherwise, 

whenever either a or b is admissible. From this equality the assertions follow as 

s I / | o otherwise. 

□ 
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From lemma|3]we easily derive a similar result regarding the block structure of the matrices obtained by commuting 
the Hamiltonians of the global gates in with the Hamiltonians of the global gates JIJ. The result can be interpreted 
as stating that, restricted to the subspace spanned of the admissible states, the commutators coincide (up to a scalar 
multiple) with the Hamiltonians of the corresponding (local) two-qubit gates acting on the pair of qubits at positions 
p and q. 

Proposition 4 For 5 G {0,1} let Us be any of the Hamiltonians -(|1<5)(0<5| - |0<5)(1<5|) and -i(\l6)(06\ + \06)(lS\). 
Let p =/= q G P, a, 6 : D — > {0, 1} such that either a or b is admissible. Then 

{a\[-iA^- r \[~iA^ r '\U^ q) ]]\b) = ~W{p,q)W(p,r)W(p,r')(a\U^ q) \b}. 

Proof. We give the proof only for Us — —\15)(0S\ + \Q8)(15\, as the calculations for the other case are essentially the 
same. Observe that Us = —Bs + B\. Hence, using also that the matrices A^ p ~ rS> and A^ p ~ r ) are self-adjoint, 



-iA^~ r \ [~iA^- r '\ U { l~ q) ]] = -[A^~ r \ [A^- r '\u'f' q) ] 



By Lemma this gives 

(a\[-iA^ r \ [-iA^~ r '\U^ q) 

whenever cither a or b is admissible. 



[A(P~ r \ [A( p - r '\B t f- q} }} - [A^- r \ [A^- r '\B { ^ q) 
[A&- r \ [A {p - r '\B [ ^ q) ]] - [A {p - r \ [A^ r '\B { v- q) ]]K 



W(p, q)W{p, r)W{p, r')(a\B ( f' q) \b) - W{p, q)W{p, r)W{p, r'){a\B { ^ q) \b) 
-~W(p,q)W(p,r)W(p,r')(a\U {p ' q ' ) \b), 



This result can be used to show that local gates on pairs of qubits in P can be efficiently simulated using global 
gates. To prove this we will need some standard facts regarding approximations of unitary operators. 

For an operator U on the Hilbert space C™ we denote by \\U\\ the operator norm of U: \\U\\ = supui_j \Ux\. Note 
that \\AB\\ < \\A\\ ■ \\B\\. If \\Ai\\, \\A 2 \\, \\Bi\\, \\B 2 \\ < 1 then we have 



\\A X A 2 - BxBiW = \\{A X - B 1 )A 2 + B X {A 2 - B 2 )\\ < \\A t - Bt\ 
By an easy induction we obtain 



\Ao 



N 



|Si||-|U 2 -S 2 || < \\A 1 -B 1 



\Ao - B 2 



\\A 1 ---A N -B 1 ---B N \\<J2\\A j -B j \\, 

i=i 

whenever At,..., An and Bt, ■ ■ ■ , Bn are sequences of unitary operators. 
Lemma 5 There is an absolute constant c > 0, such that 

\ N 

exp( ^C/- 1 )-cxp( — V- x )-exp( — U) ■ exp( —V) ) - exp([-iU, -iV\) <c-M 3 N~^ 

Vn Vn Vn Vn ) 

for any N > M 2 , where U and V are Hermitian operators on the Hilbert space C" and M = max{||[/||, ||V||, 1}. 
Proof. We use the first three terms in the expansion of exp(— tiU): 

cxp(-hU) = 1-UU- h 2 U 2 + 0(M 3 t 3 ), 
as t — > 0. The norm of the error term can be indeed upper bounded by 

oo oc 1 OO 

-{tM) j = M 3 t 3 ^ — ^{tMy < M 3 t 3 ~(tMy = M 3 t 3 e tM < e ■ M 3 t 3 



(7) 



3=3 



J- 



ft0' + 3)< 



3=0 ■ 
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if t < 1/M. Doing the same for exp(— UV), exp(UU), and exp(ttV) and collecting the terms with exponent greater 
then 2 we obtain 

exjp(-itU)- 1 ■ exp(-itV)- 1 ■ exp(-ttU) ■ exp(-itV) 
= (1 + itU - ^U 2 ){1 + itV - ^V 2 )(l - itU - yU 2 )(1 ~ itV - !jV 2 ) + 0(MH 3 ) 
= l + t 2 (VU -UV)+0(M 3 t 3 ) = 1 + t 2 [ilf,iV] + 0(M 3 t 3 ). 

On the other hand, taking just the first two term of the expansion of exp([— till, —UV]), we obtain exp([— tiU, —UV]) « 
1 + t 2 [lU, iV] , where the norm of the error term can be upper bounded by 

(2M 2 ) 2 t i e iM2t2 < 4e 4 • M 4 t 4 < 4e 4 M 3 t 3 
whenever t < 1/M. This, and the preceding formula gives 

||exp(-ztC/ _1 ) ■ exp(-i*V -1 ) ■ exp(-itU) ■ exp(-itV) - exp[-UU, exp -UV] II < c- M 3 t 3 



for t < 1/M with some constant c. Writing t = 1/yN in the latter inequality we obtain the asserted result using 
formula J7J). □ 

Now we are in a position to prove our main technical result from which theorem^ will easily follow. 

Theorem 6 Assume that for every p' =/= q' G D and p" ^g"ePU {r, r'} we have \w(p'''q")\ — w ' Then, for every 
real —1 < T < 1, for every p =/= q € P, for every S € {0, 1}, and for every < e < 1, operations which act as the two 
qubit gates 

exp(-T(|l£)((W| - |(W)(1<5|)^) , exp (-Ti(\lS){Q5\ + \Q5){l5\) {p ^ (8) 

on admissible states can be e- approximated by a product of (nw/e)°P' global gates of type £7J) and 

Proof. Let U stand for any of the global Hamiltonians i{\16){05\ - \0S) (lS\)^' q '> and (\ld)(05\ + \05) (lS\)^- q \ Let 
V = [l/W(p,r')\ll)(ll\ (p ~ r,) ,-iT/W(p,q)U]. Note that for any p, q e P, for any s € {r,r'}, \\-iT/W(p,q)U\\ < 2nw, 
1 / W(p,s)|ll)(ll| (p_s) <nw, 1/W(p,r)|ll)(ll| (p_r) < nw, \\V\\ < 4w 2 n 2 . 

By Proposition 01 we need to approximate the operation exp([— i/W(p, r)|ll)(ll|^ p ~ r - ) , — lV]). By Lemma|5]this 
can be done with error at most e/2 using a product of 0(N2) operations which are either global operations of the 

form exp(±i/(M^(p,r)v^V2)|ll)(ll|^ _r ' ) ) or operations of the form exp(±i/\/N^V), where N 2 = 0(w 12 n 12 /e 2 ). 

Furthermore by formula l|7)l. we obtain an e-approximation if we use e/(2iV2)-approximations instead of 

the operators exp(±jy / lV2^/). By Lemma we can e/(2iV2)-approximate the operators exp(±i\/~N2V) — 

exp{±[-iW(p,r'y 1 N 2 / 1/4 \ll)(ll\ {p ~ r '\-iTW(p,q)- 1 N~ 1/4: )U}) by a product of N ± global gates where N x = 
0(w 6 n 6 N* /2 e- 2 ) = 0(w 12 n 12 e~ 3 ). 

The total number of global gates used in the approximation of the "local" one is 0(N\N2) — 0(w 2A n 2A e~ 5 ). □ 

Lemma 7 The gates 

exp(-i(|L5)((W| -\06){1S\)), exp {-U{\lS){Q5\ + \0S){16\)) . (9) 
for S € {0, 1} and for real numbers —1 < t < 1 form a universal set of two-qubit gates. 
Proof. We claim that the following six 2 2 x 2 2 matrices generate SU22 as a Lie algebra over R: 

u = |ii)(oi|-|oi)(ii|, y = *|ii)<oi|+»|oi)(ii| ) 

T = |11>(10| - |10)(11|, X = »|11)(10| +*|10>(11|, 
V = 1 10) <00| - |00)(10|, Z = t|10)(00| +«|00)(10|. 

Indeed, a basis of su 2 a can be obtained as 

T, U, V X, Y, Z, 
[V,T] = |00)<11| - |11)(00|, [T,Z] = *|00)<11| +t|ll)(00|, 
[T,U] = |01)(10| - |10)<01|, [Y,T] = t|01)(10|+i|10)(01|, 
[[V,T],U] = |00)(01|- |01)(00|, [pr,z] J y]=*|oo)<oi|+*|oi)(oo| > 
[U,Y] = 2t(|ll)(ll| - |01)(01|), [T,X] = 2t(|ll)(H| - |10)(10|), [V,Z] = 2<|10)(1Q| - |00)(00|). 
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From the claim the assertion follows as the matrices U, V, Y, and Z are Hamiltonians of operations of the form J^J , 
while T and X can be obtained form U and Y, respectively, by exchanging the two qubits. □ 

We now prove our main theorem: 

Proof of Theorem ^ Consider a quantum computation (circuit) on |P| qubits. Because of Lemma [7| we may 
assume that the circuit is given as a product of £ gates of the form acting on qubit pairs in P, i.e., gates given in 
(jHJ). (The complexity in terms of other, more standard gate set is polynomially related to £.) Let e > 0. By (0, we 
obtain an e-approximation of the circuit if we use e/£-approximations of the gates. By Theorem the effect of an 
individual 2-qubit gate on admissible configurations can be approximated with error at most e/l using 0((wn£/e) k ) 
global two-qubit gates of the form and for some constant k < 24. In view of this, simulation of the entire 
circuit requires 0(£(wn£/e) k ) < 0((wn£/e) k+1 ) global gates. □ 



VI. CONCLUSION 



In the present work we have considered the computational power of a lattice composed of a two dimensional system 
(a qubit) at each site. The only gates we used were global two-qubit gates which act in a translationaly invariant 
manner on pairs of qubits. The initial state of the lattice consists of all qubits in the |0) state, except two specific 
qubits which are in the |1) state. With these ingredients we have shown that it is possible to efficiently simulate a 
quantum computer. We hope these results will stimulate further work on the computational power of lattice systems. 

Preliminary investigations suggest that one can extend the present work in several directions. First of all it should 
be possible to decrease the number of different types of global two-qubit gates which are used in the simulation. 
Secondly we have not exploited in the present work the global one-qubit gates. Preliminary work shows that they can 
be used to simplify some aspects of the simulation. 
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